Normal state of highly polarized Fermi gases: simple many-body approaches 
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We consider the problem of a single J, atom in the presence of a Fermi sea of ] atoms, in the 
vicinity of a Feshbach resonance. We calculate the chemical potential and the effective mass of 
the X atom using two simple approaches: a many-body variational wave function and a T-matrix 
approximation. These two methods lead to the same results and are in good agreement with existing 
quantum Monte-Carlo calculations performed at unitarity and, in one dimension, with the known 
exact solution. Surprisingly, our results suggest that, even at unitarity, the effect of interactions is 
fairly weak and can be accurately described using single particle-hole excitations. We also consider 
the case of unequal masses. 
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The investigation of ultracold Fermi gases with two 
unbalanced hyperfine states (which we shall denote as t 
and J.) has been through an impressive expansion last 
year. This subfield is of great interest both on practical 
and on theoretical grounds. Indeed on one hand it is re- 
lated to other fields of physics, namely superconductiv- 
ity, astrophysics and high-energy physics where similar 
situations arise On the other hand the additional 
parameter provided by the population imbalance should 
provide a tool to deepen our understanding of the BEC- 
BCS crossover in these systems and contribute to improve 
our control of many-body theory. This recent activity has 
been started by striking experimental results [2| which 
have given rise to a considerable number of theoretical 
works 0. Experiments performed on trapped systems 
have observed equal density superfluid states as well as 
partially and fully polarized regions. 

The analysis of the T = phase separation requires 
the knowledge of the properties of both the superfluid 
and the partially polarized normal phase 0,0,0. This 
has been developed in the recent work of Lobo et al Q 
where, at unitarity, the properties of the partially po- 
larized phase have been obtained by calculating through 
a quantum Monte-Carlo (MC) approach the parameters 
which characterize a single | atom immersed in a Fermi 
sea of t atoms, with density — fc|,/(67r^). In this way 
they were able to obtained a very good agreement with 
experimental results. 

Here we consider the general problem of a single J. atom 
in a completely polarized t atom Fermi sea. The t — J. in- 
teraction is characterized by an s-wave scattering length 
a whose value can be tuned via a Feshbach resonance 
from the BEC {l/kpa > 1) to the BCS {l/kpa < -1) 



limits. The Fermi gas is assumed to be ideal due to the 
suppression of higher angular momentum scattering at 
low temperatures. This problem is a much simpler one 
than the case of two equal spin populations in the BEC- 
BCS crossover, although still quite nontrivial due to the 
absence of a small parameter in the strongly interacting 
regime. It is the simplest realization of the moving impu- 
rity problem and it bears a strong similarity with other 
old, famous and notoriously difficult condensed matter 
problems, such as the Kondo problem, the X-ray singu- 
larity in metals and the mobility of ions Q and ''He 
atoms Q in '^He. This gives to the present atomic system 
a much wider significance beyond the context of polar- 
ized Fermi gases. We may hope to have a full control 
and understanding of this system, and to obtain further 
physical insight in these kind of problems. 

In this paper, in addition to considering a general scat- 
tering length a, we will extend the parameter range by 
treating the general case where the masses m| and 
are different, although they are the same in experiments 
up to now. This can be accomplished by using atoms 
belonging to different elements. The limit mi/m^ —>■ 00 
is then compared to the fully solvable problem of a fixed 
impurity in a Fermi sea. 

Here we are interested in the two physical quantities 
which have been calculated by Lobo et al [gj, namely 
the chemical potential /i| of the J, atom and its effective 
mass m* . We have addressed this problem by two differ- 
ent and complementary many-body methods. The first 
one is a natural extension of the many-body trial wave 
function used by Chevy Q to obtain the effective mass. 
The second one is a T-matrix approximation, whose ba- 
sic ingredient, ladder diagrams, appear in more elaborate 
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schemes. This approximation has already been used in a 
wide range of physical systems, including hi gh super- 
conductors M and the BEC-BCS crossover [111,113, and 
it is known to give reasonable results. The variational 
method has the advantage of providing a rigorous upper 
bound to the energy while the T-matrix can be easily ex- 
tended to include more sophisticated approximations. As 
we will see these two methods lead exactly to the same 
results and moreover, are in surprisingly good agreement 
with MC calculations [el-flst. 

The trial wave function we consider, for a system 
of total momentum p, is the following momentum eigen- 
state (we set h= \ throughout the paper): 

k>kF 

IV) = <^o|p)i |0)t + '^qklP + q - k)x 4^ cq, |0)t (1) 

q<kp 

where Ck^ and cj^.^ are annihilation and creation opera- 
tors. In the first term the f-spin free Fermi sea is in its 
ground state |0)| — nfc<fcj? "^ki bac) and the j-spin atom 
is in the plane- wave state |p)x, while in the second term 
it is in excited states corresponding to the creation of 
a particle-hole pair in the Fermi sea with momentum k 
and q respectively, the |-spin atom carrying the rest of 
the momentum. The coefficients and (/)qk are found 
by minimizing the total energy. This wave function is by 
construction suitable to reproduces the molecule in the 
BEC limit as well as the perturbative mean-field limit. 
We follow the procedure of Ref.P], in particular with 
respect to the handling of the zero range interaction po- 
tential and the corresponding regularization in terms of 
the scattering length. We obtain for the change in energy 
E due to the addition of the |-spin atom: 
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f{E,p,q) 2™ 'j,^^ eTfc + <=ipfq-k-eT<?-^ 

where ei,ik — k'^ /'^'m-ld is the kinetic energy of the f and 
J, atoms, and rrir = m|mj/(m| +mi) the reduced mass. 
For p = we have E — fi^, while the variation of E for 
small p gives the effective mass. 

These results can be obtained exactly from the knowl- 
edge of the self-energy E(p, w) of the J, atom. Indeed 
the pole of the J, Green's function Gi{p, uj) — [lu — eip + 
/i| — S(p, w)]~^, giving the dispersion relation of the J, 
quasiparticle, satisfies: 



Hp 



(3) 



Since lu is the energy measured from the chemical po- 
tential, and that physically the chemical potential corre- 
sponds to the addition of a particle with zero momentum 
p = 0, the chemical potential is given by: 



Moreover the effective mass, giving the dispersion rela- 
tion of the i quasiparticle, is obtained from the small p 
behaviour of the dispersion relation and is given by: 
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where the derivatives of the self-energy are taken for p = 
cij = (S is real in the situations we deal with below). 

The self-energy itself is obtained [l3| from the un- 
known two-particle vertex F, the only very important 
simplifying feature of our problem being that the |- 
atom Green's function is exactly the bare one, namely 
Go^{k, uj) = [u! — e^k + Mt]""'^' where /i| = kp/2m-^, since 
the single | atom does not perturb in the thermodynamic 
limit the free Fermi sea of f atoms. To obtain actual an- 
swers, we proceed to take the simplest approximation for 
the two-particle vertex, namely the T— matrix approxi- 
mation. In this approximation the | atom interacts only 
with a single t atom. This implies in particular that 
the only excited states of the | Fermi sea coming in this 
problem are the single particle-hole excitations, just as in 
Eq.([T]). In this T— matrix approximation the | and the | 
atoms scatter any number of times through the bare po- 
tential, just as in any two-body problem solved exactly 
by perturbation theory, which leads to the well-known se- 
ries of "ladder" diagrams. When there is a single f atom 
the problem reduces to the full scattering of the | and t 
atoms, and the solution can be expressed in terms of the 
scattering amplitude. When this problem is compared 
to the one in the presence of the Fermi sea of t atoms, 
the bare interaction V can be eliminated in favor of the 
scattering properties. In our approximation the vertex 
F depends only on the total momentum K and the to- 
tal energy of incomin g pa rticles. After performing the 
above steps we find [ll,]3: 
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Gi {k, n + fi^ ~ eiK-k) (e|K-k - m)] 



where 9{x) is the Heaviside function. In obtaining this 
result we have made explicit use of the fact that we will 
find < as it is obvious physically. 

In our approximation the self-energy is given by: 



2nrJ {2itYJc 



-co) (7) 



Mi = S(0,0) 



(4) 



where, as mentioned above, the fact of using Go| as the 
|-atom Green's function does not imply any approxima- 
tion. Contour C goes anticlockwise around the left-hand 
side part Recu < of the oj complex plane. We de- 
form it to enclose only the singularities of the integrant. 
Physically the singularities of T{K, fl) correspond to the 
continuous spectrum of the scattering states of the t and 
I atom, and possibly to a bound state of these atoms. 
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On the other hand Goi{k,u!) has just, for fixed fc, a pole 
at w = e|fc — /Lt|. Since we work at T = 0, the particle- 
particle continuous spectrum does not contribute as it is 
physically obvious. Indeed, from Eq.®, the correspond- 
ing singularities are at SI > (because of the Heaviside 
function) and are outside contour C. On the other hand 
we will mainly restrict ourselves to the case where there 
is no bound states between the f and the | atoms. Such a 
bound state exists clearly in the BEG limit l/kpa — > +oo 
where molecules will be present. While as we shall shortly 
show our simple approach recovers also this limiting be- 
haviour, the effect of a bound state in the intermediate 
regime will be investigated in further work. This leaves 
us only with the contribution of the pole of Got: 



S(P,W): 



dK 

(2^ 



6'(mt -eTK-p)r(ii:,cj + e|K-p-MT) (8) 



^ I dKK^ (r(K + p, + e^K - Mt)) 



where the bracket is for the angular average over the 
direction of K. 

At this stage it is interesting to consider the weak 
coupling limit a — > 0_, in which case the first term 
dominates in the right-hand side of Eq.(l6|). This gives 
T{K,Vl) = 2-Kalmr and = S(0,0) = Mi = 

kpa/ {STTrrir) — 2Tm'\a/mr. Hence we recover the ex- 
pected result for the mean-field interaction energy with 
a short-range interaction V{y) = {2'Ka/mr) We can 
view the more general result Eq.® as having a similar 
physical interpretation, but with the effective interaction 
having now a wavevector and energy dependence. 

In the general case Eq.® together with Eq.® pro- 
vides an integral equation for Y,(j),uj). We leave the ex- 
act solution of this equation for further work, whereas in 
the present paper we will proceed to a further approxi- 
mation. We will just stop at the first step of an iteration 
loop which would provide the full answer, that is we will 
set (jj) = in the expression of Gx(p, io). In this case 
we can check that Eq.Q and ([3]) are identical, provided 
E is identified with lj + ^<^. 

For the case of the J. atom chemical potential, we are 
just left with solving numerically an equation for fx^. At 
unitarity l/kpa — and for mi = this gives ui = 
—0.6066 in remarkable agreement with the QMC 0,01 
result MX -(3/5)(0.97±0.02)/i| -0.58 ± 0.01 mt. 
This surprising agreement suggests that the effect of in- 
teractions is weak even at unitarity. 

It is interesting to investigate the regime where the ra- 
tio p = Imx becomes large. In this case the expression 
of T{K, n) becomes quite simple. We find to dominant 
ordcrT{K,n) = [27r/(m^A:j.)][l/A:Fa-W{l-hr))i/2]-\ 
with r = mi/m-\^ leading to the following explicit equa- 
tion for the relation between p and 1/kpa: 
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FIG. 1: Reduced | atom chemical potential l^il as a function 
of l/kpa for various mass ratios r = mi/mi. From top to 
bottom r — 0.25 , 0.5, 1. (full thick line, the interpolating ap- 
proximation Eq.([9]) being the thick dashed-triple dotted line 
just above), and oo (lower full thick line), the exact result 
Ea. (|10|l being the dotted curve just above. The inset com- 
pares, at unitarity, the approximation Eq.([9]) (dashed line) 
with the numerical results (full line) as a function of the mass 
ratio r. 

In the case of equal masses to^ = TOj, this is plotted 
for comparison in Fig[l] and is seen to give a quite good 
agreement with the numerical value [l5|. In the weak 
coupling regime a —>■ 0- we recover the mean field result 
given above, while the asymptotic behaviour for large p 
is p = (1 + l/r)(l/fcFa)^, which is the two-body bound 
state energy. This formula can be seen as an interpolation 
between these two extremes. At unitarity it gives p — 
(1 -t- l/r)(2/37r)^/'^. For equal masses we get p « 0.71 
which is fairly near the numerical result. 

In the inset of Fig [1] we present the results of our model 
at unitarity as a function of the mass ratio r = mi/m-\. 
Again the interpolation Eq.®, i.e. p = (2/37r)2/3(l -|- 
l/r), is in quite reasonable agreement with numerical 
results. For small r the chemical potentia l ui goes to 
— cxo, as it can be seen easily from Eq.® [165. the 
other limit mi — > oo this ratio is seen to saturate. 

Assuming that, for the purpose of calculating the 
chemical potential, the thermodynamic and infinite mass 
limits commute, the problem reduces to that of an impu- 
rity interacting with a free Fermi sea, which is well-known 
in solid state physics 17[ . It can be solved exactly in the 
following way. The Fermi sea is enclosed in a large sphere 
of radius R, with i? — > cxd in the thermodynamic limit and 
the impurity at the center. Since the s-wave functions 
have to be zero at the sphere, the allowed wavevectors 
kp are given by kpR + So^kp) — pn with integer p < n 
and kpR — nn. For low energy atoms the phase shift is 
given by tanSo^k) — —ka. The energy of all the atoms 
of the Fermi sea is E — J^p fcp/2m|. The calculation is 
conveniently performed by finding the change in energy 
due to a change in the scattering length. In this limiting 
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FIG. 2: Relative effective mass m* /mi as a function of l/kpa 
for various mass ratios r = mi/mi. Same conventions as in 
FiglUfor r=0.25 , 0.5, 1. The dashed-dotted is r=4. and the 
dashed-triple dotted is r=10. The inset shows m* /mi as a 
function of r at unitarity. 

case the change in total energy of the Fermi sea is iden- 
tified with the chemical potential of the [ atom. We find 
in this way: 



\jh\ 



(1 + y^)(o+arctan2/)-t-y 



y- 
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The result is plotted in FiglT]and is seen to be in excellent 
agreement with the variational and T-matrix results. In 
particular we find p — 0.5 at unitarity, to be compared 
to the approximate value 0.465. 

Finally the relative effective mass m* /mi is displayed 
in Figl2]for various mass ratios. The first striking feature 
is that the mass enhancement is quite small around uni- 
tarity, whereas we might have expected a much stronger 
effect around resonance. Naturally when we go further 
to the BEC side m* /mi increases rapidly. Quantita- 
tively our result m* /mi = 1.17 for equal masses at uni- 
tarity is in quite reasonable agreement with the QMC 
result 1.04(3), taking into account that effective mass 
is more sensitive to approximations than energy. The 
other noticeable feature of Fig[2]is the weak dependence 
of m* /mi on the mass ratio r, as can be seen in the inset 
at unitarity. We note that, within our approximation, 
no bound state appears in the plotted range. We have 
to keep in mind that, in improved approximations, the 
location for the appearance of a bound state might be 
somewhat changed. 

To further check the reliability of our approach we have 
done the calculations in one dimension. In this case we 
can compare our results with the exact solution for equal 
masses [lq|. We find a very good agreement along the 
whole crossover, from BCS-like state to the molecular 
Tonks state, for the energy, while the mass is more sensi- 
tive to our approximations precisely when the two-body 
bound state plays a major role. 



In conclusion we have found that the physical proper- 
ties of a single | atom in the presence of a Fermi sea of 
t atoms can be described fairly accurately by the simple 
inclusion of single particle- hole excitations. Two equiv- 
alent schemes based respectively on a many-body trial 
wave function and simple T-matrix approach were de- 
veloped and give very good agreement with known MC 
results for m| = mi at unitarity. 

We are grateful to S. Giorgini, M. Yu. Kagan, X. Ley- 
ronas and S. Stringari for stimulating discussions. The 
"Laboratoire de Physique Statistique" and "Laboratoire 
Kastler Brossel" are "Laboratoires associes au Centre 
National de la Recherche Scientifique et aux Universites 
Paris 6 et Paris 7". C. L. and A. R. acknowledge sup- 
port by the Ministero dellTstruzione, deU'Universita e 
della Ricerca (MIUR). F. C. acknowledges support from 
Region He de France (IFRAF). 



[1] 
[2] 



[3] 

[4] 
[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 

[13] 
[14] 

[15] 

[16] 



R. Casalbuoni and G. Nardulh, Rev. Mod. Phys. 76, 263 
(2004). 

G. B. Partridge, W. Li, R. I. Kamar, Y. Liao and R. 
G. Hulet, Science 311, 503 (2006); M. W. Zwierlein, A. 
Schirotzek, C. H. Schunck, and W. Ketterle, Science 311, 
492 (2006)Y- Shin, M. W. Zwierlein, C. H. Schunck, 
A. Schirotzek, and W. Ketterle, Phys. Rev. Lett. 97, 
030401 (2006). C. H. Schunk, Y. Shin, A. Schirotzek, M. 
W. Zwierlein, and W. Ketterle, cond-mat/0702066 
See for example D. E. Sheehy and L. Radzihovsky, 
cond-mat/0607803 and references therein. 
A. Bulgac and M. M. Forbes, cond-mat/0606043l 

F. Chevy, Phys. Rev. A 74, 063628 (2006). 

C. Lobo, A. Recati, S. Giorgini and S. Stringari, Phys. 
Rev. Lett. 97, 200403 (2006). 

G. D. Mahan, Many-Particle Physics (Kluwer Academic, 
2000). 

N. V. Prokof'ev, Phys. Rev. Lett. 74, 2748 (1995) and 
references therein. 

F. Arias de Saavedra, J. Boronat, A. Polls, A. Fabrocini, 
Phys. Rev. B 50, 4248 (1994) and references therein. 
M. Yu. Kagan, R. Fresard, M. Capezzali, and H. Beck, 
Phys. Rev. B 57, 5995 (1998) and references therein. 
P. Fieri and G. C. Strinati, Phys. Rev. B 61, 15370 (2000) 
and references therein. 

R. Combescot, X. Leyronas and M. Yu. Kagan, Phys. 
Rev. A 73, 023618 (2006) and references therein. 
S. Pilati and S. Giorgini, unpublished. 
For an introduction to many-body techniques, see for 
example A. A. Abrikosov and L. P. Gorkov and I. E. 
Dzyaloshinkski, Methods of Quantum Field Theory in 
Statistical Physics (Dover, New York, 1963). 
Ref. [l3j] has calculated p on the a < side for equal 
masses and their results are in excellent agreement with 
ours. The agreement is not as good for the effective mass, 
but still within 10 per cent in the worst case. 
It is known that, for large mass imbalance, few-body phe- 
nomena might play a significant role in the many-body 
physics of a mixture. However here we will restrict our- 
selves to our ansatz and not take this type of correlations 



5 



into account. [18] J. B. McGuire, J. Math. Phys. 7, 123 (1966). 

[17] J. Friedel, Phil. Mag. 43 153 (1952); F. G. Fumi, Phil. 
Mag. 46 1007 (1955). 



